Photon position measure 
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The positive operator valued measure (POVM) for a photon counting array detector is derived and 
found to equal photon flux density integrated over pixel area and measurement time. Since photon 
flux density equals number density multiplied by the speed of light, this justifies theoretically the 
observation that a photon counting array provides a coarse grained measurement of photon position. 
The POVM obtained here can be written as a set of projectors onto a basis of localized states, 
consistent with the description of photon position in a recent quantum imaging proposal [M. Tsang, 
Phys. Rev. Lett. 102, 253601 (2009)]. The wave function that describes a photon counting 
experiment is the projection of the photon state vector onto this localized basis. Collapse is to the 
electromagnetic vacuum and not to a localized state, thus violating the text book rules of quantum 
mechanics but compatible with the theory of generalized observables and the nonlocalizability of an 
incoming photon. 



I. INTRODUCTION 

Although position of a massive particle is an observ- 
able in nonrelativistic quantum mechanics, photon posi- 
tion is a controversial concept. It has been argued that 
there is no photon number density, only energy density 
[l|. Solutions to the photon wave equation are electric 
and magnetic fields, but the relationship between these 
fields and the normalizable Laudau-Peierls (LP) photon 
number amplitude is nonlocal [2|, |3|. Localization of a 
converging (or diverging) photon pulse cannot be exact 
since, according to the Paley-Weiner theorem, it must 
have subexponential tails [J]. 

The most appropriate real space description of the 
photon is a subject of ongoing debate. Since photons 
are detected through interaction with matter and atomic 
dipoles are sensitive to electric field, the positive fre- 
quency part of this field is often called the photon wave 
function [5] . A LP-like wave function has been defined in 
the context of the inverse problem, spontaneous emission 
of a photon by an excited atom JQ|. A measurement of 
photon position should be described by projection onto 
a basis of localized photon states [7[, but this probabil- 
ity amplitude is not a positive frequency electromagnetic 
field. 

In spite of this theoretical controversy, a single photon 
is counted in only one pixel of an array detector and there 
is considerable interest in the measurement of photon po- 
sition in quantum optics. Position entanglement has been 
demonstrated experimentally to be a way to achieve high 
dimensional entangled states for information processing 
Q. Recently Tsang defined a photon POVM and applied 
it to a new quantum imaging method [9|. The use of a 
POVM to describe photon position was introduced first 
by Kraus llOj]. 

Use of a POVM bypasses the theoretical difficulties be- 
cause a POVM is just an operator partition of the iden- 



tity that does not require a self-adjoint position operator. 
Here a photon position POVM is obtained by lifting it 
to a larger Hilbert space that includes the atoms of a 
photodetector. We will start by considering absorption 
of photons by a two dimensional (2D) photon count- 
ing array detector. This experiment is described by a 
projection valued measure (PVM) that follows the usual 
rules for an observable in quantum mechanics. Projec- 
tion operators onto each pixel will be defined such that 
their expectation values equal the probability to create an 
electron- hole pair. We will then examine how this pho- 
ton counting experiment gives information about photon 
position. 



II. THEORY 

In the theory of generalized observables a POVM is 
an operator partition of the identity that describes the 
possible outcomes of a measurement. According to Ncu- 
mark's dilation theorem a POVM on the Hilbert space V 
can always be lifted to a PVM on a bigger Hilbert space, 
% =1^ ® A, where A is an ancillary subspace. The op- 
erators associated to a measurement onH, Fn, are then 
a set of self-adjoint projection operators satisfying the 
usual postulates of quantum mechanics. The POVM on 
V is the partial trace on the ancillary Hilbert subspace, 
A, that is 



Pn = TrA PaFu 



(1) 
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This formalism will be applied to a photon counting mea- 
surement: In Subsection A the localized states needed to 
define photon density will be discussed. In B, time de- 
pendent photodetection operators will be derived, and in 
C the special case of a thick detector that counts photons 
will be considered. Limitation of the model are discussed 
in D, and in E and F the counting operators will be ap- 
plied to one and two photon states and collapse. Normal 
incidence and the paraxial limit will be assumed so that 
cross talk between pixels can be neglected. 



A. Localized states 

A state with definite position is an equally weighted 
sum over all momenta. States of this form were used 
by Tsang to describe photon position Q. Here peri- 
odic boundary conditions over the volume V will be 
used so that the momenta, k, are discrete. Since all 
momenta are included, the two transverse polarizations 
should be defined for all wave vectors. In k-space spher- 
ical polar coordinates the definite helicity unit vectors 

ej^ ^ = ( ^ -I- ia(j)j /V^ can be defined for a — ±1. In the 

paraxial approximation the basis states are then circu- 
larly polarized (CP). The most general transverse basis 

vectors, e^l — exp (—iax) ^u i- (^' ^) ' ^^^ rotated about 
k through the Euler angle X (6*, </>) jlJl. For example, if 
-^ — —(f) the unit vectors are rotated back to the x and 
j/-axes and describe linearly polarized (LP) photons in 
the paraxial limit. 

If the operator ak,cr annihilates a photon with wave 
vector k and helicity a, the operator 



a-a (r, t) = —1= y ak,CT exp (ik • r — ikct) 
vV , 



(2) 



annihilates a photon with helicity a at position r and 
time t. The one-photon localized states are 



Ir, i,cr) 



4(r,O|0> 



(3) 



where a'l is the adjoint of aa and |0) is the electromag- 
netic vacuum state. This is an orthonormal one-photon 
basis at any fixed time, t. The number density operator 



^('''^) ^^'^li^^t)'S'^i^^t)' 



(4) 



polarization in the p-direction by a Glauber photodetec- 
tor at r will be described by the SP projection opera- 
tor |er p) {cr.pl in the localized Wannier basis satisfying 
(er,p|er',p') = SrySpy [l3|. The factor Sp^p' describes ab- 
sorption of a photon with definite polarization, as would 
be the case for a spherically symmetric hydrogen-like 
atom making a transition from a ns to a (n + l)px, 
(n + 1) py or {n + l)pz state. For normal incidence in 
the paraxial approximation we can choose z parallel to 
the direction of propagation of the beam and p = a for 
any wave vector that is actually present in the pulse. In 
this case we will only need ao- (r, f) defined by ([2]), al- 
though more generally ak,p = Y.a '^k.o- (ek, 
Gp is a mutually orthogonal triad of unit vectors 
The SP counting operators are 



e*) where 



F^P = lA®lv = l 



for zero photons. 
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(5) 



(6) 



. . (7) 

for two photons, one in pixel n and one in pixel n' ^ n, 
and so on. The sums are over polarization and all atoms 
in the n*'' pixel of the photodetector and {r G D„} are 
the positions of the atoms at which an elecron-hole (e-h) 
pair can be created. The unit operator, 1-p, acting in the 
photon Hilbert subspace V , acknowledges the fact that 
no information is obtained directly from the photons. 

These operators will be transformed to the IP and then 
to the HP. The SP Hamiltonian is 



H 



SP 



Ho +Hf'^ + Hj 



rSPt 



(8) 



B. Photodetection 

It is assumed for concreteness that the photodetector 
is made up of a 2D array of semiconductor pixels with the 
z-axis chosen parallel to its inward normal. The dipole 
moment of the individuals atoms is taken to be isotropic 
with components all equal to /i. The atoms are assumed 
to lie on a lattice so the atom positions are also dis- 
crete. The area. A, and the time, T, occupied by the 
electromagnetic pulse to be detected are taken to be fi- 
nite so that V = cTA. To obtain time dependent oper- 
ators that describe photon annihilation the Schrodinger 
picture (SP) operators will be transformed to the inter- 
action picture (IP) and then to the Hcisenberg picture 
(HP). 

Glauber defined "an ideal photon detector as a sys- 
tem of negligible size (e.g., of atomic or subatomic di- 
mensions) which has a frequency-independent photoab- 
sorption probability" |12|. Detection of a photon with 



with zero order term Hq and photon-atom interaction 
described by 

Sr = M E ^^^^'"^ (r) l^r.p) {9\ ■ (9) 

p.reV 

In ([ni, |er,p) {g\ creates an e-h pair at r with atomic dipole 
moment in the p-direction. The SP positive frequency 
electric field operator is 

^^+)'^^(r)=iE^fcexp(zkT)ak,p (10) 

k 

where, with k = Ikl, 



Ek = 



hkc 
2^' 



(11) 



In the IP and within the rotating wave approxima- 
tion Hq is unchanged but the interaction Hamiltonian 



becomes time dependent. In a semiconductor the atoms 
interact, so the Wannier states, |er,p), do not have defi- 
nite energy. The eigenvectors of the atomic part of Hq 
are the Bloch basis states, |'(/'q,p), and its eigenvalues form 
a band of energies, e^, where Cq equals to the gap energy 
plus the energy of the electron hole pair. The iV-atom 
Wannier basis is related to the Bloch basis according to 

M 



-7^'^exp{iq■r)\^Pq,p). (12) 



For a cubic lattice using periodic boundary conditions 
the wave vector q has components qi = 2'Kni/aNi for 
rii — I to Ni where a is the interatomic distance and 
NxNyNz = N is the number of atoms in the pixel. For 
the photon, the expansion pHI) gives states with definite 
energy, kc. Thus in 1^ at.p |V'q,p) (ql acquires a time 
dependence exp [i (eg — kc) t] in the IP |14j, equivalent to 

the unitary transformation O^^ = U^O^^Uq with 



Uo{t) ^exp[-iHat/h 



(13) 



We can transform from the IP to the HP by solving 
ihdU/dt = HiU iteratively to give the Dyson series |15| 



f/ (t) = 1 + f/(i) (t) + L/(2) (t) + 



where 



C/(i) (i) 



dt'HJP [t') , 



(14) 



(15) 



and it is assumed that the detector is in its ground state 
at time t^. The Hermitian conjugate has been omitted 
since photon emission by the semiconductor will be ne- 
glected. 

All photons incident on an ideal photon counting de- 
tector must be absorbed, by definition of such a detec- 
tor. Thus an J-photon counting operator must reduce a 
J-photon state to the vacuum state, |0). This is the case 
for terms of the form 



i?,f^ = c/('^)t(i)i?jpf/(j)t(i). 



(17) 



where Ai is the time required for the measurement. The 
effect of the unitary time development operators can 
be evaluated by transforming the localized basis vec- 
tors, |e 



r,pl 



to the Bloch basis vectors that are delocal- 
ized within a pixel. The Wannier and Bloch bases are 
both complete so X^r l^r,p) (er,p| can be replaced with 



Eq IV-q.p 

and 



(V-. 



q.pi 



and (112]) can be substituted in Icr 



r' ,p| 



Using assumptions equivalent to those in 
Kimble and Mandel's photodetection theory [la] in the 
impulsive detector sensitivity limit Il7| and the orthogo- 
nality of the Bloch states, a factor \L dtq exp {—iegr) ~ 
27rM (r) is obtained. With r = i" ^ 
be written as 



-t', Eq. ([HI) can then 



Fl'Iito) 



ta+At 



to 



dt Y. ^W(r,t)^(+)(r,t)|5)(g| 



P,rG-D, 



?( + ) 



(19) 
where Ep^' (r, t) is the IP positive frequency electric field 
operator, s — 2npt \^\ /h is the one-atom sensitivity to 
the absolute square of the field per unit time, and pe 
is the number of states per unit energy interval in the 
semiconductor. A Glauber ideal photodetector has a fiat 
frequency response, equivalent to assuming that s is con- 
stant, and we have followed that convention here. Ab- 
sorption of photons in different pixels is only correlated 
through the wave function, so the multi-photon counting 
operators will be time ordered integrals of factors of the 
form ([T9| . 

Inside the detector the z-component of the wave vector 
acquires an imaginary part, iak, due to absorption, and 
the speed of light becomes c/n where n is the index of 
refraction. Frequency and the component of the wave 
vector tangential to the interface is unchanged by entry 
into the detector. If a mode with free space wave vector 
k approaches the detector at an angle 9 to its normal 
it is refracted to 9d inside the detector so that sin 9d — 
si'n9/n and kz — > nkcos9d -\- iaj-. For normal incidence 
in the paraxial approximation cos 9d~\- The IP positive 
frequency electric field operator inside the photodetector 
is then 



&+'^{v,t) = iY,Eka^,^P 



X exp {ikxX + ikyy + inkz — ikct — akz) . 



(20) 



In the HP the zero photon operator remains 1, and the 
one-photon operator is 



Ki) 



|2 rto+At j-ta+At 

/ dt' / dt" 

'to Jto 



ft2 



p;r,r',r''G-DT, 



X Ul it') \g) {e,,J E^r^'^ (r') Uo (f) {e.^p) {e,,p\ 
xU^,it")El+^'^ir")\er--,p){g\Uo{t") (18) 



C. Photon counting 

The response of an ideal photon counting detector to 
photon flux is discussed by Bondurant [l3] . The detector 
must be thick enough to absorb all incident photons, so 
we can integrate over all z inside the detector. To in- 
tegrate ()19|) over z we can convert the sum over atom 
positions to an integral using ^^ ^ Pa J d'^r where pa is 



the density of atoms to give 

pto+At 

Kn' it) = / 

J to 

where A„ is the area of the n*'* pixel and 



dt I dxdywi{x,y,t)\g) {g\ (21) 

A„ 



wi {x,y,t) = spa 



p 



dz^(-)(r,t)^(+)(r,t) (22) 



and 



is the first-order coincidence rate operator in [1/ 

The essence of Bondurant's calculation can be un- 
derstood by consideration of a single mode field with 
wave vector k and definite polarization. Using (j20p and 
dnj, i?(-)i?(+) = hkca'<a/2eoV = {hkc/2eo)n where 
n = a^a/V decays inside the dielectric. Since dwi (r, t) = 
{spahkc/2eo) ndz is the absorption probability per unit 
area per unit time over the distance dz, dn = —dwi/c = 
— {spafik/2eo) ndz implies that 



ak 



spahk 
4en 



(23) 



Using Equations pT|) and ([20]) . integration of (P^ over 
a thick detector then gives wi = spa L dzE'^ E = cn^ 
that is energy density is replaced with number density, 
implying that a thick photodetector counts photons. 

More generally integration of (P^ over z with substi- 
tution of (EOl) and ^ gives [Til 



wi {x,y,t) 



E 



'kk'a 



k,p°k',p 



2eo^ 4^ , afc + afc' +in{k-k') 



(24) 



X exp \i [kx — k'x) X + i {ky — /Cy) y 
-i (A: - k') ct] . 

A photon should be detected with certainty if we sum 
over all pixels and count for a time T. To test this 
we can evaluate ^^^ Pin for t = and /S.t — T . 
Substituting ([M)) in ([T^ and summing over pixels we get 
a factor 5k^^k' Sk .k' ^k.k' inside the k and k' sum with the 
result 



n p,k 



l5) {9\ 



(25) 



Thus a photon is counted with certainty if the detector 
is big enough and it is operated over a long enough time. 
The /c-dependent factor in ([M|) can be expanded in 
powers of the relative differential photon frequency 



4 = (fc - fco) /fco 



(26) 



where fegc is the center frequency of the pulse. To first 
order in Sk 



wi (x,y,t) 



- y 

p,k,k' 



aJcp^k' 



(27) 



X exp [i (kx, — k'x) x + i [ky — k'y) y 
-i{k-k')c{t-T)] 



where r is a small delay of the order of a few optical 
cycles [ITII. This can be written as 



wi {x, y, t) = en [x, y,0,t- t) 



(28) 



where n (r, t) given by (j4|) is the number density opera- 
tor. This implies that a photon counting array measures 
photon flux density. The factor exp {—ikoct) disappears 
from ((28|) . so n describes the envelope function of the 
pulse and is unaffected by its center frequency. The z- 
integration is performed independently in each pixel and 
t > t' and t < i'are equally likely, so the two-photon 
correlation operator is 

W2 {x, y, t; x', y', t') = lyi'^^n (a;', y' , 0, t' -T)n (a;, y,Q,t- r) 

(29) 
According to ([T}, the one-photon POVM is the partial 

trace of (|2ip over the atoms, Pi.„ = Trj{ ( 'p^F^n ] ■ Fo'" 
a photodetector initially in its ground state p — p-p ® 
\g) {g\ so that the trace operation just eliminates \g) {g\ 
giving 

_ /-to+At |. 

Pi,7i{to)^ dt dxdywi{x,y,t) . (30) 

Jto J An 

If two photons are absorbed the POVM is 

P2.n,n' (to) = / dt dxdy (31) 

X / dt' dx'dy'w2 {x,y,t;x' ,y' ,t') . 

Jto J A^, 

These are coarse-grained photon position operators. The 
pixel location can be identified with the x and y posi- 
tion of the photon, while the z position of the photon is 
determined by the measurement time. 



D. Limitations of the model 

The photon counting array considered here will not 
count photons in all beams configurations. In the ex- 
treme case of propagation away from the detector, no 
photons will be counted at all. For modes not normal 
to the detector surface some photons will be counted in 
pixels adjacent to their entry pixel, and spatial resolution 
will be reduced. To go beyond the paraxial approxima- 
tion or for application to beams not normal to the de- 
tector surface, (j20p and equations derived from it would 
have to be generalized. 

A real photodetector has a frequency dependent quan- 
tum cfhcicncy less that one that should be corrected for 
when the photon number is determined from the num- 
ber of electron-hole pairs created. If it is not possible to 
determine whether more than one photon was absorbed 
in a particular pixel, the photon density should be low 
enough that multiple photon absorption in a single pixel 
is unlikely. 



It was assumed in ([6]) that an absorbed photon has 
definite polarization. This was proved to be the case for 
normal incidence with Ez neglected in [la|. In general, 

the Hermitian tensor Tr \ pEp' (r, t) E\^' (r, t) i can be 

diagonalized so "there always exist a set of three (com- 
plex) orthogonal polarization vectors such that the field 
components in these directions are statistically uncorre- 
lated" [131. If we use these polarization vectors as a basis, 
the diagonal form is exact. 

To first order in an expansion over the relative dif- 
ferential photon frequency a thick photodetector counts 
photons with a small delay time of a few optical cycles. 
The delay t has a simple physical interpretation. It is 
due to the z-integral that samples the square of the field 
over a distance of a few penetration depths, and these 
photons arrived at an earlier time. 

Second order terms in the Sk expansion that 
go as Si and 6kSk' result in uiq^q [d^a/dt^) and 



Wq {da/dt) [da/dt] in real space where wq 



cfco. 



There is also a term that goes as 5k ij^x^^l) l^o that 
gives Wq^Uq'^ (ffa/dt) [d'^a/dx^ -\- d'^a/dy'^^ in real space. 
These higher order terms depend on the rate of change of 
the pulse envelope relative to the optical frequency and 
are negligible in most experiments. 

For a single atom the factor vk gives a nonlocal 
relationship between field and number amplitude such 
that the field at r' of a photon at r is proportional to 
|r — r'|~ ' @;ll3l- This no nlocality is almost completely 
eliminated by integration over z, and is neglected here. 



E. Photon state and wave function 



As an illustration ([28| and ([29l) will be applied to 
paraxial one and two photon states at normal incidence. 
For a one-photon pulse, if the time required for a mea- 
surement is less than the pulse length, detection of either 
1 or photons is possible. The probability to count one 
photon between times fo and to + ^i is Pi = X^n ^i," 
given by (I30|) while the zero-photon probability is 1 — Pi . 
If this intermediate photon number is read out, the wave 
function collapses to a or 1 photon state [l^. To 
avoid the complications of intermediate collapse we will 
only consider measurement times greater than the pulse 
length. 

A HP one-photon pure state is of the form 



1^1 



E 

k,(T 



Ck,, 



(32) 



The photon wave function will be defined here as the 
projection of j-^i) onto the one-photon state ([3]) with po- 
larization a localized at r at the measurement time i, 
that is 



V'o-(r,i) = (r,i,cr|V^i) 

Eexp (— ik • r -|- ikct) 

k 



V 



(34) 



At each point in space ip„ has two components corre- 
sponding to the two transverse photon polarizations. For 
a CP paraxial beam the unit vectors at = reduce to 
ej_i = (x -I- zcry) /\/2 and, together with z, can be iden- 
tified with the Gp basis. For LP we can take x = — f/* so 



that e^-*^ 



X and e 



y at 



0. In either case 



(■01 \wi{x,y,t)\il)i) =c^|-0<^(a;,y,O,i)|^ (35) 
<j 

is photon flux density. 

A two-photon HP state vector is of the form 

IV'2) = ^ Ck,cr;k',<T' |k, (T;k',CT') . (36) 

k,k',(j,cr' 

If we define the symmetric two-photon wave function 

V^,,,. (r, t; r', t') = ^ (0| S, (r, t) Z,, (r', t') 1^2) (37) 



'2\ 



= pr-^ / . (Ck,(T;k',(T' + Ck\cr';k,<T) 

X exp [i (-k • r - k' • r' + feci + fc'ci')] , 

the two-photon coincidence rate is 

(V'2|w2(a;,y,i;a;',y',i')IV'2) 

= c2^|Va,.'(a:,y,0,i;x',y',0,Ol'- (38) 

For example, if the photons are produced by sponta- 
neous parametric down conversion (SPDC) it was found 
that their positions are highly correlated 8]. For type 
II SPDC the polarization of the photons is orthogonal, 
while for type I the polarizations are the same. Both 
these effects can be described by ([55]) . In the photon 
counting experiment considered here polarization was not 
measured, and should be summed over. Knowledge of 
photon polarization comes from knowledge of the beam 
properties rather than from the photodetection process. 



where X^k a- I'^k,^! — 1- The expectation value of the 
photon flux density operator is 



(■01 |wi(a;,y,t)|0i) 

= c ^ (-01 ^\ (x, y, t) I 0) (0 |a^ (x, y, t) \ -01 



(33) 



F. Collapse 

We will start by considering a pure state, |0i) \g), that 
initially contains one-photon and assume that an e-h pair 
has been detected in pixel n. The extended system of 
photons plus detector atoms satisfles the usual rules of 



quantum mechanics. In the IP this state vector evolves 
with time to * (t)) = U^^^ (t) \^pl) \g) and collapses to 
the final state ^(^) — F(^^ |* (t)) after a measurement. 
Since its final and initial states are the same, F/^ — Ff^ 



is given by (J6l) so that 

E K.)(er,.|C/W(t)|^i)|ff). (39) 



i*r/) 



oc 



(T=±l,reD„ 



If Tpa (x, y, 0, t) is uniform over the area of a pixel and po- 
larization isn't measured the normalized collapsed state 
vector simplifies to 



\Kf) 



1 



(T,re-D„ 



(40) 



where TV is the number of atoms in the n pixel. If J 
photons are detected and the wave function is approxi- 
mately uniform over a pixel the system collapses to the 
state 



^ 



)=n 



1 



E 



(41) 



<Ji,rieD„ 



III. DISCUSSION 

The one-photon fiux density operator (l^5|) for he- 
licity a can be written as c |r, ijCr) (r, i,(T| evaluated 
at z = 0. The Hermitian position operator r^ — 
J (Prv |r, t, a) (r, i, a\ with eigenvectors |r, t, a) and eigen- 
values r and a can be defined such that the usual 
rules of quantum mechanics are formally satisfied. For 
two photons the projection operators can be written as 
|r,i,(T; r',i', cr') (r, ijCr; r',i', cr'l, consistent with the ID 
multi-photon POVMs defined by Tsang 9]. There are 
in principle many ways to construct a POVM, since any 
partition of the identity is allowed. However, an exactly 
localized position eigenvector is an equally weighted sum 
over all wave vectors, and this form arose naturally here. 
But there is an important deviation from the usual rules; 
collapse is not to this position eigenvector. 

By the usual rules of quantum mechanics, if position 
of a photon with polarization a is measured at time t the 
system should collapse to the instantaneously localized 
state, |r, i, a) . Since all wave vectors are equally likely in 
a localized state such as ^ , both incoming and outgoing 
states are required [7|. In a photon counting experiment 
there is no outgoing wave and collapse to an exactly local- 
ized state is not possible. In place of the usual rules, here 
we have used the theory of generalized observables that 



separates the calculation of probabilities from collapse. 
With this formalism we found that a photon counting 
array detector measures coarse grained photon position, 
while collapse is to the zero photon vacuum state. Some 
theorists would insist we say "absorption position" was 
measured, but the term "photon position" is widely used. 
Definite helicity basis states were selected because all 
wave vectors can be included and helicity is Lorentz in- 
variant. States with a definite polarization direction in 
real space exclude longitudinal wave vectors and thus can 
not be localized. The definite helicity unit vectors have 
definite total angular momentum parallel to some z-axis 



11[. Selecting the z-axis is analogous to selecting the 
quantization axis for electron spin that is arbitrary in 
the absence of symmetry breaking, but can be fixed by a 
magnetic field. For the photon, symmetry can be broken 
by the beam axis. 

In nonrelativistic quantum mechanics the real space 
wave function is both {v\^j) and a solution to the par- 
ticle's wave equation. For photons these are not the 
same thing, the former being a scalar while the latter 
is a vector whose k*'' term is weighted as y/k. Here 
the photon wave function is defined in the former sense, 
that is as the projection of the photon state vector 
onto a basis of the localized states. The photon wave 
function is i/'cr(i",i) = (r,i, o'|'i/'i) and the photon flux 
density measured in a photon counting experiment is 
c^^ IV'ct {x, y, 0, t)| . The quantum electrodynamic state 
vector is used to describe the photon state and a wave 
function that satisfies Maxwell's equations is not re- 
quired. This choice of wave function is consistent with 
Chan et al [61 , with Tsang [9| , and with the discussion of 
position eigenvectors in |7|. 

In summary, a photon position POVM for a photon 
counting experiment has been derived by performing a 
partial trace over the ancillary Hilbert space of the pho- 
todetector atoms. It was found that the photon counting 
probability equals the photon fiux density integrated over 
the area of the pixel and the time taken for a measure- 
ment. This is, of course, as expected, since a photon 
counting array is designed to measure photon position. 
If the photon wave function is defined as the projection 
of its state vector onto the localized states, photon fiux 
density is its absolute square multiplied by the speed of 
light. The photons are destroyed in this measurement, 
but the distribution of positions and times at which they 
are detected is a measure of their probability density in 
the incoming pulse. 
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